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The high pt > 3 GeV azimuthal asymmetry, V2(pt), in non-central nuclear collisions at RHIC is 
shown to be a sensitive measure of the initial parton density distribution of the produced quark-gluon 
plasma. A generalization of the Gyulassy-Levai- Vitev (GLV) non-abelian energy loss formalism 
including Bjorken 1+1D expansion as well as important kinematic constraints is used. 
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Introduction. In order to interpret data on nuclear 
collisions from recent Relativistic Heavy Ion Collider 
(RHIC) experiments [ [l], ^, ||, it is obviously necessary 
to have knowledge of the initial conditions. Currently, 
there is an order of magnitude uncertainty in the ini- 
tial produced gluon density, p g (T ) ~ 10 — 100/fm 3 , in 
central Au + Au at s/s — 130 AGeV since widely dif- 
ferent models [ g, §] are consistent [ § with PHOBOS 
data [ [J . We note that recent PHENIX data [ |J appear 
to be inconsistent with at least one class (final state [ 
El) of gluon saturation models. It is essential, however, 
to check this with other observables as well. High p^ 
observables are ideally suited for this task because they 
provide a measure [ Q of the total energy loss, AE, of 
fast partons, resulting from medium induced non-abelian 
radiation along their path [ ffl, H|- For intermediate jet en- 
ergies (E < 20 GeV), the predicted [§, |o) gluon energy 
loss in a static plasma of density p g and thickness, L is 
approximately AE GLV ~ E(L/6 fm) 2 / 9 9 /(10/fm 3 ). The 
approximate linear dependence of AE on p g is the key 
that enables high observables to convey information 
about the initial conditions. However, AE also depends 
non-linearly on the geometry, L, of the plasma and there- 
fore differential observables which have well controlled 
geometric dependences are also highly desirable. 

A new way to probe AE in variable geometries was 
recently proposed in Ref . [ [ll| . The idea is to exploit the 
spatial azimuthal asymmetry of non-central nuclear col- 
lisions. The dependence of AE on the path length L(<p) 
naturally results in a pattern of azimuthal asymmetry of 
high px hadrons which can be measured via the differen- 
tial elliptic flow parameter (second Fourier coefficient), 
V2(pt) [ 0. In this letter, we predict V2(pt > 2 GeV) 
for two models of initial conditions [ ^| which differ by 
an order of magnitude. We first generalize the finite en- 
ergy GLV theory [ H| to take into account the expansion 
(neglected in [ [0j [n|) of the produced gluon-dominated 
plasma while retaining kinematic constraints important 
for intermediate jet energies. Another novel element of 
the analysis is a discussion of the interplay between the 
azimuthally asymmetric soft (hydrodynamic [ [l^]) and 
hard (quenched jet) components of the final hadron dis- 
tributions. We show that the combined pattern of jet 
quenching in the single inclusive spectra and the differen- 
tial elliptic flow at high p^ provide complementary tools 



that can determine the density as well as the spatial dis- 
tribution of the quark-gluon plasma created at RHIC. 

Hadron transverse momentum distributions. It is use- 
ful to decompose the nuclear geometry dependence of in- 
variant hadron distributions produced in A + B ^ h + X 
at impact parameter b into a phenomenological "soft" 
and perturbative QCD (pQCD) calculable "hard" com- 
ponents as 

dN AB (b) = N part {b) dN s (b) + T AB {b) da h (b) , (1) 

where N part (h) is the number of nucleon participants, 
and T AB {h) = / d 2 r T A {r)T B {v - b) is the Glauber pro- 
file density per unity area in terms of nuclear thickness 
functions, T A (r) = J dz p A (r,z). The computable low- 
est order pQCD differential cross section for inclusive 
p + p — ► h + X production is given by 

i pp /> 
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where x a ,x b are the initial parton momentum fractions, 
z c = Ph/Pc is the final hadron momentum fraction, 
fa/p(%a,Qa) are the parton distribution functions and 
Dh/c(z c ,Qc) is the fragmentation function for c — > h. 
The UA1 data on pp hadron production with p^ > 1 GeV 
can be well reproduced with the above formula us ing 
Q 2 = p|/2, K = 2 and Martin-Roberts-Sterling [ |f 
(MRSD-') structure functions. 

In nuclear collisions jet quenching can modify the hard 
cross section by changing the kinematic variables of the 
effective fragmentation function. We follow Ref. [ and 
include this effect by replacing the vacuum fragmentation 
function in Eq. (0) by an effective quenched one 



z cD' h/c {z c , Ql) = z' c D h / c (z' c , Q 2 C ,) + N g z g D h/g (z g ,Q 2 g ) 
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where z' c , z g are the rescaled momentum fractions. The 
first term is the fragmentation function of the jet c after 
losing energy AE c (p c , (f>) due to medium induced gluon 
radiation. The second term is the feedback due to the 
fragmentation of the N g {p c ,<j)) radiated gluons. The 
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modified fragmentation function satisfies the sum rule 
J dz c z c D' h/c (z c ,Q 2 c ) = 1. 

Energy loss in a longitudinally expanding plasma. The 
GLV reaction operator formalism [ |9j expands the radia- 
tive energy loss formally in powers of the mean number, 
X, of interactions that the jet of energy E suffers along 
its path of propagation through dense matter. For a jet 
produced at point xq, at time To, in an expanding and 
possibly azimuthally asymmetric gluon plasma of den- 
sity p(x, r), the opacity in direction v(4>) is 



x(<t>) 



dr o-(t)p(x + v{4>){t - t ), r) 



(4) 



Note that the gluon-gluon elastic cross section, cr(r) = 
97rag/2/igyj(r), and the density may vary along the 
path. For a finite jet energy, E, the approximate upper 
kinematic bound of medium induced momentum trans- 
fers is |q(r) max | » y/3p(T)E and Me//( T ) = M 2 ( T )(1 + 
A <2 (' r )/ t l 2 (' r )max)- The explicit closed form expression 
for the n th order opacity expansion of the gluon radi- 
ation double differential distribution for a static medium 
is given in Ref. [ |). Fortunately, the opacity expansion 
converges very rapidly due to the formation time physics, 
and the first order term was found to give the dom- 
inant contribution. Higher order corrections decrease 
rapidly with energy. All numerical results in this let- 
ter include 2 nd and 3 rd order correction factors com- 
puted in the static plasma limit [ H]. We also include 
finite kinematic bounds on the transverse momentum, 
kmax = mm [^E 2 x 2 , AE 2 x(l — x)] and k 2 lin = /i 2 , for glu- 
ons with light-cone momentum fraction x. Finite kine- 
matics reduces energy loss at intermediate jet energies [ 
!§, [l(| as compared to the asymptotic formalism [ gf. 

The dominant (generalized) first order radiation inten- 
sity distribution [ |9(] that holds also for expanding plas- 
mas is given by {z = r) 
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In order to compare to previous asymptotic results [ 
for expanding plasmas, consider a density of the form 



p(z)=p Q {^ ) a 0{L-z) 



(6) 



where a = corresponds to a static uniform medium of 
thickness L, while a = 1 to a more realistic Bjorken 1+1D 
expansion of the plasma (transverse to the jet propaga- 
tion axis) . Analytic expressions can be obtained only for 
asymptotic jet energies when the kinematic boundaries 
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can be ignored [ |8|] . If we set q; 
gleet the running a s and change variables k — q — > k, 
u = k 2 /p 2 (z) and w — q 2 / ' p 2 (z) 1 then Eq. (E) reduces to 



dx 



E 



2C R a s 



dza(z)p(z) f(Z(x,z)) 



(7) 



where Z(x, z) = (z — zo)p 2 (z)/2xE and 
du 



f(x,z) 



u(l 



[1 — cos ( uZ(x, z))] 



^ + ^-lo g (Z)+0(Z 2 ) 
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For a target of thickness L, the small Z(x, z) limit ap- 
plies as long as i > i c = Lp 2 (L)/2E. In that do- 
main dl/dx oc 1/x. For x -C x c , f{Z) w log Z and 
dl/dx cx log l/x is integrable to x = 0. 

By integrating over x, the total energy loss is 



C R a s 



dZ ; ; ; [Z — Zq) log 



X(z) 



p{z) ' 



(9) 



which is an approximately linearly weighed 
line integral over the local transport coefficient 
(p 2 (z)/X(z)) log E/n(z) = 9iTa 2 p(z)v/2. For a uniform 
and expanding plasma as in @ 

„2-q - 
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Here v = logE/p and we used that p? (L) L a / X(L) is 
a constant independent of L for this type of expansion 
and took the Zq — * limit. We therefore recover the 
asymptotic Baicr-Dokshitzer-Mueller-Schiff (BDMS) and 
Zakharov (Z) energy loss for both static and expanding 
media [ p|. We note that for Bjorken expansion, the 
asymptotic energy loss can be expressed in terms of the 
initial gluon rapidity density as 



AE a=1 (L) 



1 dN9 \ L E 
ttR 2 dy ) ^ p 
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If we vary L = R oc A 1 / 3 by varying the nuclear size, 
then nonlinearity in L arises because dN 9 /dy oc A 1+s . 
For HIJING initial conditions [ [§ 5 = 1/3, while in the 
EKRT saturation model [ § 6 « 0. 

Implications of nuclear geometry. For nucleus-nucleus 
collisions the co-moving plasma produced in an A + B 
reaction at impact parameter b at formation time t = zq 
has a transverse coordinate distribution given by 

1 da iCt 

Pg (v, z = 0,t = z o ) = — T A (r) T B (r - b) , (12) 

z dy 

where da ict /dy is the pQCD mini-jet cross section in pp 
collisions at a given y/s. Note that taking into account 
also the 2D transverse expansion causes the density to de- 
crease somewhat faster than Eq. (||). However, we found 
numerically that transverse expansion can be ignored in 
the first approximation. 
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FIG. 1. The GLV fractional energy loss in Bjorken expand- 
ing gluon plasma with [ [| § dN 3 jdy ~ 200, 1000. 
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FIG. 2. The modulation function i?(b, </>) is plotted vs. <f> 
for several impact parameters. and Wood-Saxon vs cylinder 
geometries. 



In the linear f(Z) « irZ/2 (||) and Bjorken approx- 
imations, the total energy loss is proportional to the 
line integral (Q,|9|) along the jet trajectory r(z, (f>) — 
r + v((f))(z — zq), averaged over the distribution of the 
jet production points 



F(b,cj)) = / (Tr 



j2 T A (r)T B (r-b) 



TAB(b) Jza 
T A (r(z,<f>))T B (T(z,<l>)-b) , 



(13) 



where T^r), T B (r — b) and T^sih) depend on the geom- 
etry. In particular, for a sharp uniform cylinder of radius 
i? e ff one readily gets Ta(t) = (A/nR 2 s )9(R c ff — |r|) and 



?ab(0) = A 2 /nR^ s . We can therefore define the effec- 
tive radius of the sharp cylinder equivalent to a diffuse 
Wood-Saxon geometry via 



F(0, 



) Wood — Saxon 



F(0,, 



/Sharp cylinder 



(14) 



For Au + Au collisions and a — 1, Eq. ( |l4| ) gives 
RcS ~ 6 fm. Eq. (^|) can then be integrated numeri- 
cally to give AE(0)/ allowing a s to run and including 
kinematical bounds. Fig. 1 illustrates the fractional en- 
ergy loss for gluon jets at b — for a broad range of 
initial gluon densities [ |[ [|. 

For a non- vanishing impact parameter b and jet direc- 
tion v {<()), we calculate the energy loss as 



AE(b, 
E 



F(b, t 
F(0,< 



AE(0) 
E 



(15) 



where the modulation function i?(b, (f>) captures in the 
linearized approximation the b and <f> dependence of the 
jet energy loss. Fig. 2 shows the R(b, 4>) modulation 
factor plotted against the azimuthal angle cf> for impact 
parameters b = 2, 6, 10 fm. Note that i?(b, 4>) reflects 
not only the dimensions of the characteristic "almond" 
cross section shape of the interaction volume but also the 
rapidly decreasing initial plasma density as a function of 
the impact parameter. 

Phenomenological soft "hydrodynamic" component. In 
order to compare to the new STAR data [ |^| at pr < 
2 GeV, we must also take into account the soft non- 
perturbative component that cannot be computed with 
the eikonal jet quenching formalism above. In [ [ll| this 
was simply modeled by an azimuthally symmetric expo- 
nential form. However, in non-central A+B reactions the 
low pt hadrons are also expected to exhibit azimuthal 
asymmetry caused by hydrodynamic like flow effects [ 
|l2| . We therefore model the low soft component here 
with the following ansatz: 

dNJb) dn s e - pT / To , , . . ,.. 



where we take To ~ 0.25 GeV and in corp orate the hy- 
drodynamic elliptic flow predicted in [ [12[ and found to 
grow monotonically with as 

v 2s (pt) w tanh(p T /(10 ± 2 GeV)) . 



(17) 



It is important to emphasize that hydrodynamic flow was 
found [ |l2| to be less sensitive to the initial conditions 
than the high px jet quenching reported here. 

With the inclusion of this non-perturbative soft com- 
ponent, it follows from Eq. ([T]) that the effective differ- 
ential flow is 
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FIG. 3. Sensitivity of central-collision inclusive hadron dis- 
tributions to initial conditions and energy loss in the two 
component hydrodynamic + quenched jet model. 
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FIG. 4. The interpolation of V2(pt) between the soft hydro- 
dynamic [ [l^] and hard pQCD regimes is shown for 6 = 7 
fm. Solid (dashed) curves correspond to cylindrical (Wood- 
Saxon) geometries. 
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V2s(PT)dN s + V2h{pr)dN h 
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(18) 



This interpolates between the hydrodynamic and the 
pQCD regimes because at high pt, dNh ^> dN s . For our 
numerical estimates the low pt interpolation is achieved 
by multiplying the pQCD curves with a switch function 
[1 + tanh(3(p T - 1-5 GeV))]/2. 

Conclusions. Fig. 3 shows the inclusive charged parti- 
cle transverse momentum distribution in central Au + Au 
collisions with three models of initial gluon densities 
dN3/dy = 1000,500,200. We see that jet quenching can 
be disentangled from the soft hydrodynamic component 
only for transverse momenta pt > 4 GeV. In that high 
Pt region there is an approximately constant suppression 
relative to the unquenched (HIJING) distribution due to 



the approximately linear energy dependence [ |Toj of the 
GLV energy loss [ ^] . The suppression increases system- 
atically with increasing initial plasma density and thus 
provides an important constraint on the maximum initial 
parton densities produced in b = collisions. 

Fig. 4 shows the predicted pattern of high pt 
anisotropy. Note the difference between sharp cylinder 
and diffuse Wood-Saxon geometries at b = 7 fm, the 
characteristic impact parameter of minimum bias events. 
While the central (b = 0) inclusive quenching is insensi- 
tive to the density profile (due to Eq. (fl4|)), non-central 
events clearly exhibit large sensitivity to the actual dis- 
tribution. We checked numerically that transverse ex- 
pansion with v± = 0.5c can be ignored since it reduces 
the jet quenching effects in Figs. 3,4 by <20% at high px- 

We conclude that V2(pt > 2 GeV,b) provides essen- 
tial complementary information about the geometry and 
impact parameter dependence of the initial conditions in 
A + A. In particular, the rate at which the v-i coefficient 
decreases at high is an indicator of the diffuseness of 
that geometry. 
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